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Abstract. We present a new chaotic system of three coupled ordinary differential equations, limited to 
quadratic nonlinear terms. A wide variety of dynamical regimes are reported. For some parameters, chaotic 
reversals of the amplitudes are produced by crisis-induced intermittency, following a mechanism different 
from what is generally observed in similar deterministic models. Despite its simplicity, this system therefore 
generates a rich dynamics, able to model more complex physical systems. In particular, a comparison with 
reversals of the magnetic field of the Earth shows a surprisingly good agreement, and highlights the 
relevance of deterministic chaos to describe geomagnetic field dynamics. 
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■ 1 Introduction 

Deterministic chaotic models have been largely studied in 
, the past few decades, both on the numerical and theoret- 

■ ical sides [I], (2, [3]. In addition to their interest to non- 
] linear mathematics, these models often play an important 

■ role in the comprehension of physical systems. Indeed, in 
many situations, they offer a simple modelisation of exper- 

, imental and natural objects exhibiting complex dynamics. 

■ Particular interest has been devoted to models involving 
, only three equations, the minimum number required to 

■ generate chaos in a continuous dynamical system with 

■ time independent coefficients. Similarly, ordinary differ- 
, ential equations of these types of models are generally re- 

■ stricted to the lowest order terms. This approach ensures 
I that numerical integration yields chaotic dynamics rele- 

■ vant to describe the behaviors of complex systems, mean- 
I while a theoretical study of the model is still possible. 

In this article, we focus on the dynamics of reversals, 
i.e. the periodic or chaotic changes of sign of the ampli- 
tudes in dynamical systems. Reversals between two sym- 
metric states are observed in several natural or experi- 
mental systems, among which we can mention thermal 
convection, magnetohydrodynamics dynamos, or turbu- 
lent atmospheric flows. Although these systems generally 
involve turbulent states and a very large number of de- 
grees of freedom, several sets of deterministic equations 
have been proposed to model reversals, the Lorenz model 
providing the most famous example [1]. As typical exam- 
ples of chaotic three-dimensional systems exhibiting re- 
versals, we can also mention the Rossler attractor [5], the 
Rikitake equations [6] and the Sprott models [7]. Such 
deterministic three-mode equations often provide relevant 
toy models, but generally fail to describe the details of the 



dynamics of the modeled system. In the case of reversals 
occurring in fluids mechanics, in which turbulent fluctu- 
ations are important, explanations in terms of stochastic 
processes are thus usually preferred to nonlinear chaotic 
dynamics. Here, we propose a new model for reversals, 
given by the following system of three coupled ordinary 
differential equations: 

Q = HQ- VD (1) 
D = -vD + VQ (2) 
V ^ r -V + QD (3) 

We will see that the type of dynamics generated by these 
equations is relatively different from previous nonlinear 
chaotic models, and is sometimes more similar to stochas- 
tic systems. It has been recently proposed as a model for 
reversals of a turbulent magnetic dynamo obtained in nu- 
merical simulations [5]. In the last section of this article, 
we will discuss a physical interpretation of this model in 
the framework of reversals of the Earth magnetic field. 

In this system of equations, the nonlinear terms are re- 
stricted to quadratic order, and the parameters u and F 
are real. Relative signs of the quadratic terms are such that 
the solution remains bounded. The divergence of this vec- 
tor field is equal to fx—v—l, so it is dissipative for ji > 1+v. 
This very simple dynamical system presents similarities 
with several three-modes models previously discussed in 
the literature. In the context of geodynamo reversals, Rik- 
itake [5] proposed a model close to ((T]-[21), except that the 
linear damping of the mode V was discarded. Moreover, 
in Rikitake model, D and Q were linearly stable. Nozieres 
[9j derived a similar model, also related to geomagnetic 
reversals, but obtained from a truncation of MHD equa- 
tions, the physical interpretation of this model being a 
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Fig. 1: Numerical integration of the model for ^ = 0.119, ly = 0.1 and F = 0.9. Poincare section used for iterated map 
is shown in red. Top panel: trajectories in the phase space (D, Q, V){lcit) and in the cut {D, Q) (right). Bottom panel: 
Time evolution of amplitudes Z?(left), Q (middle) and l^(right). 



system of two magnetic components coupled to one veloc- 
ity mode. Therefore, in Nozicres model, linear terms rep- 
resent ohmic decays and the modes are linearly damped, 
strongly modifying the dynamics. Finally, the dynamical 
system (fD— TSj) with = is relevant for several hydro- 
dynamical problems, and has been extensively studied in 
a series of articles by Proctor and Hughes [TO] . Again, we 
will see that a non-zero F is crucial for the dynamics de- 
scribed in the present paper. Finally, note that the system 
possesses the symmetry (D, Q, V) — > {—D, —Q, V), which 
is crucial for the dynamics of reversals reported in section 
3. 



2 Stability and fixed points 

In this section, we first describe the fixed points of the 
system. A first trivial fixed point O is related to vanishing 
amplitudes of D and Q: 



Qo^Q 



Dn^O 



Vo=F 



(4) 



Linearization around this fixed point give the following 
eigenvalues: 
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(5) 
(6) 



The trivial eigenvalue A3 indicates a stable direction, inde- 
pendently of the values of /i, v and F, and traces back to 
the linear damping of the amplitude V. When fi + i^ > 2F, 
all eigenvalues are real. When ^ + v< 2F, eigenvalues Ai 
and A2 become complex conjugates and the system under- 
goes Hopf bifurcation for fi > ly. Note that this system ex- 
hibits a co-dimension two bifurcation point for fi = v = F, 
i.e. the system is close to both a stationary and a Hopf bi- 
furcation. We will see that the behavior of our dynamical 
system strongly depends on the stability of O. 

Apart from this trivial point, there are four non-trivial 
fixed points Mzp and Nzp . Two of these solutions only exist 
when -v//!^ — F > 0: 



UI± 



F>M± 



Fr- 



(7) 



(8) 
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The other fixed points are 



(9) 



{N± 



D 



N± 



(10) 



(11) 
(12) 



Linearization of the system around the fixed points 
M± gives the following characteristic relation: 



,(j^-m) 



A-' + A^(l + i.-^) + A r 



Similarly, for N± we have: 



A-' + A'(l + j^-^)-A r 



(13) 

A^v + Ar^/JIu = 0. 



(14) 

In the following, we describe some of the chaotic be- 
haviors generated by the system of equations pi— T^)) . by 
exploring the parameter space while F is kept fixed (in 
practice F = 0.9). We consider only situations where the 
linearization around non-trivial fixed points gives one real 
eigenvalue, the two others being complex conjugates. The 
real eigenvalue is always negative, ensuring that the fixed 
points have a stable direction. The real part of the com- 
plex conjugate eigenvalues will depend on the value of /.i 
and v. 

Among the wide variety of dynamical regimes observed 
in this simple system, a particularly interesting chaotic 
behavior can be observed for some parameters (see for in- 
stance the chaotic attractor shown in figure[Tl obtained for 
= 0.119, V = 0.1 and F = 0.9). For these parameter val- 
ues, the system explores the phase space around one of the 
fixed points N for a given time, and chaotically switches 
from one solution to its opposite, leading to chaotic re- 
versals of amplitudes D and Q {V does not change sign). 
Although the evolution during a given polarity is strongly 
chaotic, reversals exhibit very robust characteristics. In 
order to understand the mechanisms leading to such dy- 
namics, in particular how the system becomes chaotic, it is 
useful to keep parameters F and v to these values [F = 0.9 
and V = 0.1), while only ^ is varied. In addition, the di- 
mension of the problem can be reduced by means of a 
Poincare map in the phase space [5]. As shown in figure 
[TJ Poincare map is defined as the surface orthogonal to 
the (D, Q) plane which contains the line joining the two 
fixed points Nj^ and (this line provides a graduated 
axis with O the origin). Our first return map is then de- 
fined by the intersection Xn+i between trajectories of the 
system (projected in the (Z?, Q) plane) and this line, as a 



function of the previous intersection Xn . Figure [3] shows a 
bifurcation diagram of x„ as ^ is decreased from 0.15 to 
0.1. 
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Fig. 2: Top: Bifurcation diagram of the return map Xn as 
a function of fi. When ^ is decreased from 0.15 to 0.1, the 
transition to chaos is achieved through a period doubling 
cascade. Bottom: Evolution of the Lyapunov exponent A 
of the system when fi is decreased. The successive period 
doubling events lead to a chaotic state, characterized by 
an increase of A. 



When fj, > 0.2, two stationary solutions are obtained, 
corresponding to the fixed points N+ and N- . For ^ < 0.2, 
these fixed points loose their stability and a stable periodic 
orbit is created around each fixed points, corresponding to 
stable solutions in the iterated map. Below fi = 0.138, a 
series of pitchfork bifurcations is observed in the return 
map of figure [21-top, due to the successive period doubling 
bifurcations of the periodic orbit surrounding A^-|_ . Because 
of the symmetry of the system, a similar transition occurs 
around N- . 

We thus see that our simple dynamical system exhibits 
a classical route to chaos, by period doubling bifurcations. 
This scenario is observed in several discrete applications, 
like the Henon map or the logistic map. Note that the tran- 
sition to chaos is very sharp, with a Lyapunov exponent 
A (characterizing the rate of separation of infinitesimally 
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Fig. 3: Numerical integration of the model for /i = 0.12, 
V = 0.1 and F — 0.9. These parameters correspond to 
a fully chaotic state, with trajectories fluctuating around 
one of the fixed points. Top: discrete map. Bottom: time 
evolution of D. Initial conditions are chosen in the basin 
of attraction of iV+ . 



close trajectories) which rapidly tends to 0.1 at /i w 0.13 
(see figure [2]-bottom) . 

Figure [3] shows trajectories for ^ = 0.12, for which 
a chaotic attractor is obtained. On the iterated map, the 
system randomly explores the interval 1.4 < Xn < 2.3 (fig- 
ure 131 top) , defining a typical size of the chaotic attrac- 
tor, bounded in the phase space (D, Q, V). It is important 
to note that symmetries of our dynamical system impose 
that a similar attractor is obtained in the vicinity of the 
opposite fixed point N^. However, for /i = 0.12, these 
two symmetric attractors are disconnected, and there are 
no reversals (figure [31 bottom) . This is a noticeable dif- 
ference with other deterministic models for reversals. For 
instance, in the Lorenz model, a period doubling route to 
chaos can also be observed [TT], but in that case, the ini- 
tial periodic orbit is a cycle surrounding the fixed points 
of the system, and therefore already involving reversals 
of the amplitudes. In the next section, we describe how a 
further decrease of n yields chaotic reversals of D and Q, 
very different from reversals obtained in Lorenz or Riki- 
take models. 
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Fig. 4: Dynamics of the model for fj, ~ 0.119, v = 0.1 
and r — 0.9. Discrete application (top) is very similar to 
the one obtained for ^ — 0.12, except for a small region 
connected to the opposite fixed point. Bottom: evolution 
of the system in the phase space (D, Q). 



For /I = 0.12, the system chaotically fluctuates around 
a fixed point, but without reversing. From this state, a 
slight decrease of fi yields a transition to a dynamics of 
reversals between two symmetric states with well defined 
values. Figure |H shows the evolution of the model for 
fx = 0.119, in such a reversing regime. Note that the be- 
havior close to one of the fixed points is similar to the one 
generated for = 0.2, i.e without reversals. However, for 
H = 0.119, the two attractors are now connected, allowing 
chaotic reversals between the two opposite states. 

Changing one of the parameter of the system (here 
fj.) yields a discontinuous change of the structure of the 
attractor. This typical behavior of dissipative systems is 
known as a crisis phenomenon. The crisis behavior of dy- 
namical systems has been discovered and extensively stud- 
ied by Grebogi, Ott and York [T^], [13], mainly in the con- 
text of discrete maps. Following these authors, a crisis is 
defined as "a collision between a chaotic attractor and a 
coexisting unstable fixed point or periodic orbif^ . 
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Typically, in a crisis scenario, modifications of the struc- 
ture of the attractor are of three kinds: in a first case, 
the chaotic attractor is suddenly destroyed at the critical 
point. In a second case, the size of the attractor undergoes 
an abrupt increase. In a last case, two attractors or more 
collide and a larger attractor is created. In each of these 
three cases, the structural change is related to a collision 
with an unstable attractor (chaotic or not), and it is pos- 
sible to characterize the phenomena with a typical time r. 
This r measures the finite time during which the system 
behaves as before the crisis, and then rapidly leaves this 
region of chaotic transient. 

In our dynamical system, the apparition of chaotic re- 
versals seems to belong to the third class, i.e. the merging 
of two distinct attractors. Let us define as the criti- 
cal value of fi for which reversals appear. When /x > /ic, 
the two attractors close to each of the fixed points pos- 
sess their own basin of attraction, separated by a bound- 
ary basin. Depending on the initial conditions, trajecto- 
ries tend to one attractor or the other, and are trapped 
in the basin. Because of the symmetry of the system, 
when jjL = ^c, the two attractors simultaneously collide the 
boundary basin separating them, which allows a connec- 
tion between both attractors, through some unknown un- 
stable periodic orbits of the attractors. As a consequence, 
for fj, = fj,c, the system spends a long time exploring chaot- 
ically the vicinity of one of the fixed point before it sud- 
denly leaves this region and fluctuates around the oppo- 
site solution. The two attractors being symmetrical, the 
process then repeats and chaotic reversals are generated. 
We can thus characterize the dynamics by a time r which 
is naturally defined as the average duration between two 
successive reversals, hereafter referred as Waiting Time 
Duration (WTD). 

In general, Grebogi et al have shown that r follows a 
power law 

T^{^^~ Aie)-^ (15) 

Using numerical simulations, we determined numeri- 
cally this law. For each value of /i, t is obtained by aver- 
aging all the WTDs obtained in a time series containing 
one thousand reversals. 

Figure [5]-top shows the averaged waiting time dura- 
tion T between two successive reversals as a function of fi. 
The curve is diverging for /Xc ~ 0.11928. In practice, this 
critical value has been obtained by noting that no rever- 
sals occurred after an integration time t > 1.10^, whereas 
reversals were observed for fi = 0.119275. 

Figure [5]-bottom shows how the waiting time durations 
evolve as the critical point is approached. There is little 
doubt that this curve follows a power law. However, even 
if our numerical simulations suggest a critical exponent 
7 ~ 1/2, it is difficult to conclude on a precise exponent, 
since a small error on the determination of fic significantly 
modifies the measured power law. 

In some idealized cases, Grebogi et al. [13] have shown 
that it is possible to obtain the critical exponent 7 theo- 
retically and to express it as a function of the eigenvalues 
of the periodic orbits involved in the crisis. In the case of 
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Fig. 5: Evolution of the characteristic waiting time du- 
ration between two successive reversals when the critical 
point of the crisis is approached. Top: r as a function of 
fi, suggesting = 0.11928. Inset shows 1/t(/x). Bottom: 
Evolution in logarithmic scale as a function oi fi — iic- Nu- 
merical simulations suggest an exponent close to 7 ~ 1/2. 

discrete maps: 

7= i + (ln|ai|)/ln|a2| (16) 

where ai et 012 are respectively eigenvalues of the unstable 
and stable manifolds of the two periodic orbits considered. 
In principle, it should be possible to use this formula to 
obtain 7 theoretically from the Poincarc map. Unfortu- 
nately, the method involves calculating numerically the 
stable and unstable manifolds of the orbits mediating the 
crisis, leading to an error bar on the estimate of 7 too large 
to be compared to the exponent computed from the wait- 
ing time durations. Although the value of 7 seems very 
close to the theoretical lower bound, our numerical results 
are not accurate enough to conclude. 

This relation between the exponent and the eigenval- 
ues of the periodic orbits relies on the assumption that 
the crisis is of the heteroclinic tangent type, which means 
that it is induced when the unstable manifold of a pe- 
riodic orbit (A) becomes tangent to the stable manifold 
of a second unstable periodic orbit (B) (located outside 
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Fig. 6: Basin of attraction of chaotic attractors projected in the Poincare section. White and black colors distinguish 
the basins of attraction of each attractors. Attractors are shown in red. Top: before the crisis, for /z = 0.1194 (left) and 
fi = 0.1193 (right). Bottom: For = 0.11928 « /Zc, we observe a collision between attractors at the basin boundaries, 
illustrating the boundary crisis. 



the basin of attraction before the crisis) . One particularity 
of our dynamical system is that the orbit {B) destroying 
the first attractor belongs to the basin of attraction of 
the second attractor. Because of the symmetries, identical 
connections are achieved with both attractors, allowing 
the system to chaotically come back to the attractor it 
previously leaved. 

From the Poincare section defined in figure [1] it is pos- 
sible to visualize the collision leading to the crisis, and 
figure [6] shows the basin of attraction of the system in 
the Poincare section (X, Y) . Theses basins are numeri- 
cally obtained by using a grid of (500 x 500) initial condi- 
tions {Xo,Yo): the point {Xo,Yo) is colored in black (re- 
spectively in white) if the computed trajectory starting 
from this initial condition ultimately tends to the attrac- 
tor surrounding N- (respectively Figure IHl^a) shows 
the system before the crisis, for /i = 0.1194 > fic. In red 
are points corresponding to the intersection of trajecto- 
ries of the attractor with the Poincare section. Note that 
the two attractors arc distinct. Black and white regions 
clearly indicate symmetric basins of attraction, and are 
such that the two attractors are connected only to their 
own basin. Figure El^b) corresponds to a slightly smaller 



value, fi = 0.1193 > fic, showing that a change of fi toward 
fic yields a change of shape of the basins of attraction and 
an increase of their size. Finally, for /i = 0.11928 w /ic, 
each of the chaotic attractor collides with the boundary 
of the basin of attraction of the symmetric attractor. A 
look at the trajectories in the phase space indicates that 
the orbit involved in this boundary crisis seems to be an 
unstable periodic orbit encircling the trivial point O. Note 
indeed that this kind of reversals are only observed when 
O is an unstable focus. It is interesting to note that when 
these reversals appear, the largest Lyapunov exponent is 
not modified, since reversals are due to the fusion of two 
identical chaotic attractors nearly unchanged by the crisis 
event. In this regime, we thus obtain a Lyapunov expo- 
nent A ~ 0.1. 

As shown above, for F = 0.9 and = 0.1, a decrease of 
fi leads to four well defined successive behaviors: a station- 
ary state, periodic then chaotic oscillations around this 
non-trivial fixed point, and finally, chaotic reversals be- 
tween the vicinity of the two fixed points when fi < fXc- 
For further decrease of fi, the averaged waiting time be- 
tween reversals strongly decreases and the system finally 



Christophe Gissinger: A new deterministic model for chaotic reversals 



7 



reaches a regime of chaotic oscillations. Figure [7] she 
the evolution of the system as a function of time and 
phase space, for fj, = 0.101, v = 0.1 and F = 0.9. Whei 
finally reaches the value of j/, we note that the system ^ 
comes linearly stable at the origin O, since the real pe -S 
of Ai and A2 vanish. However, this bifurcation appeart 
be strongly subcritical: for initial conditions chosen ins < 
the basins of attraction of the non-trivial fixed points, 
trajectories stay inside this basin of attraction. The s 
tcm stabilizes on the trivial fixed point O indepcnder 
of initial conditions only for fj, < 0.063. 
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Fig. 7: Top: phase space for fi = 0.101, v = 0.1 and F = 
0.9. Reversals are not chaotic anymore for these values. 
A periodic connection is achieved between the two fixed 
points. Bottom: Time evolution of the three amplitudes 
D, Q, V. 



Inside this interesting regime of bistability, a series of 
transitions is observed in the dynamics of our system as 
fi is decreased. As it can be seen in figure [51 for 0.065 < 
^ < 0.1, the system explores the vicinity of one of the 
fixed point during a given time and suddenly makes an 
oscillation around the trivial point O. Although it stays 
on the basin of attraction of the same attractor, the tra- 
jectory seems to briefiy explore the opposite state. During 
the time spent close to the fixed point, the behavior can 
be slightly chaotic. We thus observe an intermittency be- 
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Fig. 8: Evolution of the amplitude D for different values 
of /i. Top, left: ^ = 0.09. Top, right: ^ = 0.08. Bottom: 
fi — 0.065. Note the asymmetric behavior of the dynamics 
in the latter case. However, the symmetric attractor can 
be obtained with different initial conditions. 



tween a chaotic state around a fixed point and an unstable 
periodic orbit: as time evolves, the system switches chaot- 
ically from the vicinity of the fixed point to the unstable 
periodic orbit encircling the origin. As expected, the mean 
time between these switches increase when decreases. 
For fi = 0.065, surprising asymmetric reversals can be ob- 
served, despite the symmetry of equations. 

The set of parameters chosen so far ensures the exis- 
tence of only two of the four non-trivial fixed points. We 
have seen that this is already sufficient to generate several 
different dynamical regimes, mainly related to reversals. 
One interest of this model is its ability to describe very 
different regimes, without any needs of nonlinear terms 
of higher degrees. This can be illustrated by choosing pa- 
rameters such that ^/JIl' — F > 0, meaning that there are 
now four non-trivial fixed points in addition to O. In this 
case, the two additional solutions induce new dynamics. 
Figure [9] shows the evolution of trajectories in the phase 
space {D, Q) (top) and as a function of time (bottom), for 
H = iy^ 0.35 and F = 0.1. 

Finally, others exotic behaviors can be obtained. For 
instance, figure [TO] shows the evolution of the system for 
yU = 0.3, u = 0.2 and _r = 5. In this regime, an intermit- 
tency between two distinct behaviors is observed. During 
a first phase, trajectories slowly escape from the unstable 
periodic orbit surrounding the trivial point O. When the 
amplitude of the oscillation is sufficiently large, it reaches 
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Time [arb. units] 

Fig. 9: Reversals obtained ioi fi = v = 0.35 and F = 0.1, 
when the four non-trivial fixed points are involved. Top: 
evolution in the phase space {D, Q). Bottom: evolution of 
the amplitude _D as a function of time. 



the basin of attraction of one of the two fixed points, and 
after some transient wandering around this fixed point, 
the trajectory is re-injected toward the center of the ini- 
tial periodic orbit, and the process repeat. Note that here, 
chaos is essentially due the apparent random choice of the 
fixed point involved at the end of the oscillations around 
O. 



Compared to previous similar models, this simple dy- 
namical system, involving only 3 quadratically coupled 
equations, generates several different dynamical behav- 
iors, depending on the values of the free parameters 
V and r. It would be interesting to see at which point 
the different characteristics described here are robust. For 
instance, one could add multiplicative or additive noise to 
equations ([D- TS)) , or take into account terms of higher de- 
grees to couple the three amplitudes. The present form of 
the model (without noise and limited to quadratic terms) 
presents of course the advantage to remain extremely sim- 
ple, and yet generating complex dynamical regimes. 
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Fig. 10: Evolution of the system for /j, — 0.3, v — 0.2 and 
r ^ 5. Top: evolution in the phase space {D,Q). Bottom: 
evolution of the amplitude D as a function of time. 

4 Physical interpretation and comparison 
with geodynamo reversals 

The perspective of a description of the reversals of the 
Earth's magnetic field in the framework of a low dimen- 
sional model is very tempting. Indeed, a low dimensional 
description provides a simple vision of the phenomena, 
and allows to make interesting predictions concerning the 
shape and the statistical properties of the geomagnetic re- 
versals. In fluids dynamos, the small ratio between kine- 
matic viscosity and magnetic resistivity (small magnetic 
Prandtl number), together with the proximity to the dy- 
namo threshold, strongly reduce the degrees of freedom 
of the system, and low-dimensional dynamics can be ex- 
pected. For this reason, this approach has been followed 
by several studies in the past [S] , [B] , [H] , [H] , [S] ■ Unfor- 
tunately, despite their chaotic behaviors, low dimensional 
deterministic models generally yield too simple dynam- 
ics, and fail to realistically reproduce observational data. 
For instance, Lorenz and Rikitake models exhibit typi- 
cal growing oscillations between reversals, characteristics 
of their low-dimensionality, but absent from geomagnetic 
measurements. 

As shown in the previous sections, reversals of our sim- 
ple model are generated by crisis-induced intcrmittency, a 
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Fig. 11: Comparison between time series obtained from paleomagnetic data (top, right) (SINT2000 data [17]) and from 
the model {D, Q, V) with /i = 0.119, ly = 0.1 and F = 0.1 (top, left). A typical time series obtained from the Lorenz 
model is shown at the bottom. 



mechanism different from the one creating chaotic rever- 
sals in classical models. We will show in this section that 
the model ([Tl- T^l) leads to a better agreement with gcody- 
namo reversals than Lorenz and Rikitake models. More- 
over, using symmetry arguments, it is possible to give a 
simple physical interpretation to the three amplitudes. 

The magnetic field of the Earth is mainly dominated 
by its dipolar component. Figure [TTTtop-right) shows the 
evolution of the Earth dipole moment during the past two 
millions years fl7| . Although its amplitude is not known 
during reversals, it has been suggested that the quadrupo- 
lar mode could also play an important role during rever- 
sals [IH], [H]. In the model pH- T^I) . amplitudes D and Q 
can be regarded as dipolar and quadrupolar components 
of the magnetic field. Note that the system is invariant 
with respect to the symmetry {D —D, Q — Q), in 
agreement with the symmetry {B ^ ~B) of the magnetic 
field B in the equations of magnetohydrodynamics. An- 
other important symmetry in the Earth core is the mir- 
ror symmetry with respect to the equatorial plane. Un- 
der this symmetry, the dipole is unchanged whereas the 
quadrupole is transformed in its opposite. This implies 



that the dipole and quadrupole can not be linearly coupled 
if the flow in the Earth core is mirror symmetric. By not- 
ing that the amplitude V is coupling D and Q in equations 
((T]- T5|) . we can therefore interpret V &s a. velocity mode 
(externally forced by P 7^ 0) which breaks the mirror sym- 
metry of the Earth core flow. The term QD in equation |3] 
thus stands for the Lorcntz force (quadratic in the mag- 
netic field B). The system of equations pH- TSI) therefore 
models the interaction between dipolar and quadrupolar 
magnetic modes coupled by a mirror-antisymmetric veloc- 
ity mode, and can be directly compared to geomagnetic 
data. 

This interpretation of the model is in fact very close to 
a recent model proposed by Petrelis and Fauve, also rely- 
ing on a dipole-quadrupole interaction |20| , and aiming to 
model the dynamics observed in the VKS (Von Karman 
Sodium) experiment, a laboratory turbulent fluid dynamo 
PTj . The VKS experiment exhibits chaotic and periodic 
reversals of the polarity of the magnetic field, resulting 
from a surprisingly low-dimensional dynamics [22] . and 
very similar to the predictions of the Petrelis-Fauve model. 
In this model however, chaos is generated by the addition 
of a stochastic noise in the equations (modeling the tur- 
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bulcnt fluctuations), and the velocity mode V is not ex- 
plicitly taken into account. Nevertheless, it is interesting 
to note that most of the regimes of the VKS experiment 
can also be reproduced with our simple model, despite its 
deterministic dynamics [53]. This underlines the fact that 
at small Pm, the characteristics of the turbulent flow and 
the details of the configuration (such as boundary condi- 
tions) can be ignored as far as the dynamical regimes of 
the magnetic field are considered. 

Kono has shown in 1987 that the Rikitake model was 
unable to reproduce the Probability Density Function (PDF) 
of waiting times between reversals, or the PDF of the am- 
plitude of the dipolar field [SI] . In general, distributions of 
the averaged waiting time between reversals are assumed 
to be poissonian for the geomagnetic field . It is however 
difficult to conclude, since only a small number of reversals 
are available [55]. The Rikitake model clearly fails to de- 
scribe a poissonian process: depending on the parameters, 
the model predicts either a cutoff forbidding long events, 
or the apparition of too many long events. We have com- 
puted the distribution of the waiting times for fi ~ 0.119, 
v = 0.1 and F = 0.9, and figure [T51 shows a comparison to 
the paleomagnetic data, obtained from the CK95 database 
[55] . One can see that the very short times are character- 
ized by a quantification. This obviously comes from the 
minimal time it takes to the system to follow at least one 
periodic orbit around a given fixed point before reversing. 
This behavior is typical of models with low-dimensional 
chaos. For larger waiting times however, a poissonian law 
is obtained, with a distribution 

P(t) = -^—e-<^ (17) 

< T > 

Our model is thus able to reproduce the poissonian behav- 
ior observed for the Earth magnetic field, despite the small 
number of degrees of freedom. Because the quantification 
at small time is a strong property of deterministic chaos, 
it would be interesting to search for such a behavior in 
geomagnetic observations. Note however that this type of 
behavior is easily destroyed in the presence of noise, which 
cannot be avoided in the Earth core due to the turbulent 
fluctuations of the flow. 



From a general point of view, we see that this new 
model allows to reproduce some of the main characteris- 
tics of the geodynamo, like the poissonian process or the 
distribution of the intensity of the field. As it appears in 
figure (TT] this resemblance between the two systems is 
also apparent from a simple examination of the time se- 
ries: during a stable polarity period, both the model and 
the Earth magnetic field show some chaotic oscillations, 
with sometimes bigger events called excursions, during 
which the dipole seems to reverse but finally comes back 
to its initial value. Moreover, this model reproduces the 
tendency of the Earth dipolar field to systematically over- 
shoot its quasi-stationary value just after each reversal, 
and also the time scale separation between the reversals 
and the stable polarity states. This model also makes some 
interesting predictions concerning the geodynamo. For in- 
stance, when the dipole D goes to zero during a reversal, 
this corresponds to a strong increase of the quadrupo- 
lar component Q in the model (see figure [H]) . Note that 
because of the deterministic nature of the model, the re- 
versals always exhibit these features. This predicts that 
during a reversal, the magnetic field of the Earth would 
change shape (by transferring the magnetic energy from 
the dipolar to the quadrupolar component), rather than 
simply vanishing. 

In addition, this model relies on the existence of a ve- 
locity mode V breaking the mirror symmetry of the Earth 
core fiow. Several studies seem to confirm this role of an 
equatorially antisymmetric velocity mode in the geody- 
namo reversals. First, it has recently been observed that 
the ends of superchrons (large period of time without geo- 
magnetic reversals) correspond to major fiood basalt erup- 
tions due to large thermal plumes ascending through the 
mantle [57] . In the framework of our model, these plumes 
correspond to an increase in the symmetry breaking of the 
fiow {r is increased). In agreement with this observation, 
it has been shown in global numerical simulations that 
taking an heterogeneous heat flux (rather than homoge- 
neous) at the core-mantle boundary of the Earth strongly 
influences the frequency of magnetic field reversals |28j . 
Again, this corresponds in the model to change the forc- 
ing of the antisymmetric velocity mode V. 



Similarly, there is a strong disagreement between ob- 
servations and classical models on the distribution of the 
intensity of the dipole. The geodynamo exhibits two sym- 
metric maxima around a minimum at zero, but Rikitake 
and Lorenz models statistically spend too much time close 
to zero, giving a different distribution. A study of the dis- 
tribution of the amplitude D in our model during reversals 
shows a bimodal distribution with a very low probability 
at the origin, in agreement with geomagnetic observations. 
This success of the model is mainly due the the type of 
chaos generated in the model, which produces a time scale 
separation between the duration of a reversal (related to 
some unstable periodic orbit) and the waiting time dura- 
tions (controlled by the distance from the crisis bifurcation 
point). 



5 Conclusion 

In this paper, we have presented a new deterministic model 
for reversals. In this model, chaotic reversals are generated 
through a boundary crisis between two symmetric chaotic 
attractors. Compared to previous deterministic models, 
this mechanism produces a relatively complex behavior 
in a given polarity, while reversals exhibit a very robust 
shape. In addition, we have shown that several very dif- 
ferent dynamical regimes can be obtained, depending on 
the values of the free parameters of the system. 

This model has been compared to reversals of the Earth 
magnetic field, and a good agreement has been obtained. 
Although geomagnetic reversals are certainly more com- 
plicated than a simple interaction between three modes. 
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Fig. 12: Left: Distribution of the waiting time durations (WTD) for the model. We observe large oscillations at small 
time corresponding to the quantification of the unstable periodic orbits, and a poissonian behavior at large time. 
The red curve shows an exponential fit. Right: Probability distribution of waiting times between reversals for the 
geomagnetic field, showing an approximative poissonian law (from CK95 database [26)). 




Fig. 13: Left: distribution of probability of amplitudes D. Right: bimodal probability distribution of the intensity of 
the Earth dipolar field. 



it is interesting to obtain a simple way to describe this 
complex phenomenon. The strong similarities between our 
model and paleomagnetic observations suggest an inter- 
pretation of geomagnetic reversals as resulting from an in- 
teraction between the dipolar and the quadrupolar compo- 
nents of the magnetic field. In addition, our model predicts 
that a velocity mode could play the role of the amplitude 
V, by coupling the two magnetic components dipole and 
quadrupole. It would be interesting to study the differ- 
ent possibilities for such a fiow, corresponding to a break- 
ing of the mirror symmetry of the large scale fiow in the 
Earth core. Because of the fast rotation of the Earth, 
the geostrophic equilibrium makes difficult a spontaneous 
breaking of this symmetry. However, the presence of the 
core, decoupling the north and the south hemispheres in- 
side the tangent cylinder, makes possible the existence of 
large scale antisymmetric velocity modes: for instance, an 
heterogeneous heat flux at the core-mantle boundary may 
force the large scale velocity mode V. 



Several physical systems, different from the dynamo, 
display chaotic reversals of an observable between two 
symmetric states. In general, the interpretation in terms 
of a few coupled equations generating deterministic chaos 
is easily ruled out. Indeed, most of the models predict a 
behavior far too simple between reversals, like the growing 
periodic oscillations of the Lorenz or Rikitakc equations. 
It is possible to attenuate this behavior by adding noise 
to the system, taking into account higher order terms, or 
increasing the dimension of the model. The new model 
presented in this article shows that it is in fact possible to 
obtain a more complex behavior between reversals, closer 
to stochastic models predictions, with only three coupled 
equations limited to quadratic orders. This is an impor- 
tant observation, because it opens the perspective to un- 
derstand the dynamics of complex physical systems using 
very simple mechanisms, mainly relying on symmetry ar- 
guments. 
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Fig. 14: Evolution of the three components of the model 
during a reversal, for = 0.119, v = 0.1 and F = 0.9. 
Note the rapid time scale related to the reversal event, 
compared to the rest of the dynamics. When the dipole 
D goes to zero during the reversal, the magnetic energy 
is transferred to the quadrupolar component of the mag- 
netic field. After the reversal, the dipole systematically 
overshoots its mean value. 
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